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Introduction. ~ The Heisenberg chain of n spins ^ with 


the Hamiltonian (known as the XXX model) 


n—1 

H = + 1 ), 

21=0 


( 1 ) 


where Ux = (cr^, cr^, tr^) are Pauli operators and periodic 
boundaries are assumed (?„ = Sq, is arguably the sim¬ 
plest nontrivial interacting quantum many-body model. 
The spectrum and eigenstates of H can be formulated 
in terms of the famous Bethe ansatz [l[, which gave 
birth to the theory of quantum integrable systems [2 Q . 
Eq. (IT|) has been originally prraosed as the model of 
(anti) ferromagnetism in solids and is indeed a ver; 
good description of the modern spin-chain materials 
It may also be considered as a fundamental paradigm of 
quantum statistical mechanics which is being used for de¬ 
veloping theoretical mechanisms of non-equilibrium dy¬ 
namics and thermalization or relaxation to the general¬ 
ized Gibbs ensemble (GGE) [^-Q. 

The relaxation dynamics based on quantum quenches 
[Ml gave firm evidence that the full set of (~ n) local 
conserved operators, the existence of which is granted 
for a quantum integrable system, is incomplete, in the 
sense that it cannot describe the steady state completely 
through a GGE. Similarly, a numerical experiment count¬ 
ing the number of linearly independent time-averaged lo¬ 
cal operators HI indicated that the set of local conserved 
charges should be incomplete and numerical approxima¬ 
tions of new quasilocal operators have been put forward. 

In this Letter we explicitly construct new families of 
non-local but quasilocal operators by composition of a 
transfer matrix (TM) - in the sense of algebraic Bethe 
ansatz, but for higher integer or half-integer auxiliary 
spins s > i - and its derivative, at a special combina¬ 
tion of spectral parameters, which in the thermodynamic 
limit (TL) becomes equivalent to a logarithmic derivative 
of TM. Furthermore, we prove quasilocality (in full rigour 
for a finite set of auxiliary spins s) as well as linear in¬ 
dependence of these new operator families w.r.t. local 
conserved charges. Generally, we identify quasilocality 
with the condition of factorizability of the largest eigen¬ 
value of an auxiliary TM that enters in the computation 
of the norm of the conserved operator, i.e. a product of 


higher-spin TM and its derivative. As we facilitate finite¬ 
dimensional unitary representations of quantum or Lie 
symmetries, the new quasilocal operators are always spin- 
reversal symmetric unlike in alternative recent construc¬ 
tions in the XXZ chain (l3 - [l7j| which only work at spe¬ 
cial, commensurate values of anisotropy. This features 
promise that our technique shall be applicable for gener¬ 
ating new quasilocal charges in other generic integrable 
models with Lie or quantum group symmetries. Being 
able to construct a complete, or as large as possible set 
of independent quasilocal conserved charges is crucial for 
any application in quantum statistical mechanics, besides 
constructing GGE, e.g. in linear response theory at hnite 
temperatures, studies of quantum ergodicity and many- 
body localization, etc. Quasilocal conservation laws are 
also closely related to boundary driven/dissipative quan¬ 
tum chains IJ, [181 . 

Transfer matrices and conserved operators.- Let Vg, 
s G denote a 2s-|-l dimensional spin-s module, Vg = 
= lsp{|m),m = —s, —s-I-1,..., s}, Isp denoting 
a linear span of a set of vectors, carrying the unitary 
irreducible representation of SU{2) with generators 


s^\m) = m|m), s^|m) = ■\/(s -I- 1 ± m){s =F m)\m ± 1). 

( 2 ) 

The physical Hilbert space is an n—fold tensor product of 
fundamental representations Up = V® 2 i cr^ = 2s^, 
= ^{o'^ ± icT^) = s^. Fixing arbitrary s € and 
considering another, auxiliary Hilbert space Ha = Vg, we 
dehne Lax matrices as operators over Hp G Ha 

La;,a(A) = A1 -|- tT^S^ + -|- Cr,j, s/ = A1 -|- fTx ' Sa, (3) 


where A S C is the spectral parameter. Throughout the 
Letter, operators acting nontrivially over the auxiliary 
Hilbert space are written in bold, or double strike font if 
acting over multiple (tensor product of) auxiliary spaces. 
As a simple consequence of Yang-Baxter equation, the 
(physical) TMs Ts{X) £ End('Hp) 

Tg(A) = traLo,a(A)Li,a(A) • • • L„_i,a(A), (4) 


where s is the auxiliary spin, form a commuting family 


[rg(A),Tg/A')]=0, Vs, s', A, A'. (5) 
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The fundamenal TM Ti(A) is generating the complete 
set of local conserved Hermitian operators 

n—1 

Qt = -idt" logTi (i + it)|t=o = E 0 qk ), 

(6) 

k > 2, with <52 = H, where qk G End(V® 2 ) ^~Point 
operator density, and 5 is a cyclic lattice shift map over 
End(Hp) defined by 5 (ct“) = a“^d(x+i.n)- 

Locality and quasilocality.- The 4"-dimensional space 
of physical operators End('Hp) is turned into a Hilbert 
space by defining a Hilbert-Schmidt (HS) inner product 
iA,B) := (A^B) w.r.t. the infinite-temperature state 
(H) := 2“”trH. Let {H} := A — (A)l denote the trace¬ 
less part of an operator. One of physically most impor¬ 
tant features of the local conservation laws Qk is the ex- 
tensivity of the HS norm ||{(3fe}||HS := {{Qk},{Qk}) = 
(2~tr {q^qk) — | 2 “'"'tr n oc n. We define (equiva¬ 
lently to [l^) a general traceless, translationally invari¬ 
ant operator A = iS(H) G End('Hp) as quasilocal if two 
conditions are met: (i) ||^||hs (^0 1 °" 

cally supported fc-site operator b = 6^0 l 2 n-fc the overlap 
( 6 , A) is asymptotically, as n —>■ oo, independent of n. 

One should stress that quasilocality only makes sense 
in TL n —>■ oo as it is the property of an infinite se¬ 
quence of operators labelled by n, rather than opera¬ 
tors for any fixed size n. More intuitively, a quasilo¬ 
cal operator Q can be thought of as a convergent sum 
of local operators Q = where includes 

only terms supported on r contiguous sites and the sum 
IIQIIhs = 1 II<5^''^IIhs is ramdly (typically exponen¬ 

tially) converging. Usually |14l4l6| . quasilocality can be 
detected by inspecting the leading eigenvalue of a certain 
auxiliary transfer matrices, whose r-th power yields the 
partial norm of the r—site terms The effect of 

quasilocal conserved operators to statistical mechanics is 
arguably as important as that of local operators. In par¬ 
ticular, quasilocal charges can be understood as those 
conserved operators of one-dimensional systems which 
can influence equilibrated (steady-state) values, say af¬ 
ter a quantum quench, of strictly local observables. Our 
central result is the following 

Theorem: Traceless operators Xs{t), s G G K., 

defined over the physical Hilbert space TLp as 

X,{t) = [T.(t)]-" {r,(-i + it)T;(i + it)} , (7) 

Tsit) =-t'^ - {s + , ( 8 ) 

where r}(A) = d\Ts{X), are quasilocal for all s,t and 
linearly independent from {Qk', k > 2} for s > i. 

The fact that Xs{t) are exactly conserved and 
[Xs{t),Xs/{t')] = [Xs{t),Qk] = 0 follows directly from 


The form of our ansatz 0 is inspired from an ob¬ 
servation (see Eq. m or, e.g., Ref. [l9|) that at s = 5 , 
TM becomes in TL n —>■ 00 a unitary operator 


I °° t’^ \ 

ri(i-Hit) ~ exp gfc+ij , (9) 

and hence 0 can be associated with a logarithmic 
derivative via T'J'(A) = T’J(A) = (—l)”Ts(—A) where the 
last equality is due to spin reversal symmetry 
s=*= —>• —s^. 

Proof of quasilocality.- First, we write a matrix prod¬ 
uct form of a general product of a pair of TMs 


r.(M)T,(A) = tr,„a. n 


( 10 ) 


x=0 \a£j 


where the operators L“(/i, A), a G J := {0,x, y,z} act 
over a pair of auxiliary spaces "Hai 0 "Has = Vs 0 Vs 

L°(//,A) = A//1 + Sai • ^ 2 , (11) 

l-'s (M; a) — i X Sa 2 “h ASaj^ H” 

Identity component can be written with the Casimir op¬ 
erator C = (sai -H SaJ^ as L° = ^A1 -H i(C - - sfj, 

hence its spectrum reads r|(^. A) = —s{s + l)+fiX, 
j = 0,1,..., 2s. Placing the spectral parameters along 
one of the two lines 

V^ = {{yLt,Xt)',tGR}cC\ 

+ xt:=±x+it, (13) 

we define the restricted auxiliary operators as L^“(t) := 
L“ (fif, Xf). The dominating eigenvalue of Hermitian op¬ 
erator L+°(f) = L 7 °(t) is Ts(t) = T°{p.f,Xf), Eq. ®, 
corresponding to the singlet eigenstate 

S 

IV'o) = ( 2 s + l)-i /2 ^ 0 I - m), (14) 


with a finite gap to the subleading eigenvalue 'r'(t), 5 = 
log |Ts(t)/T'(f)| > 0, for any t. The condition (sa^ + 
Sa 2 )|' 0 o) = 0 and the SU{2) algebra Sa^ x Sa^ = isa^ 
imply the following useful identities 

EWI'^o) = 0, (lpo\t+it) = -2(V'o|Sai, 
(ipolKit) = 0, L“(t)|V'o) = -2saJV'o)- (15) 

We proceed by constructing a TM over a 4-spin auxiliary 
space Ha = ^t=l ’^ak, Hai ,2 = Vs, Hi,,, ,, = Vs' 

Ts,s'(/i,A,Ai',A')= El:(m,A)0L“,(^',A'), (16) 

a^J 

which helps us to compute a general inner product of 
the form Ks,s'{t,t') := {Xs{t),Xs>{t')). The Hilbert- 
Schmidt kernel (HSK) then immediately follows after dif¬ 
ferentiating traces of powers of suitable TMs 
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= [Tsit)Ts'{t')] "'d^-d^+ (tr ,Xt - tr [L°(/^t , )]” tr [L°,(/x+, . (17) 


As a consequence of boundary condition given by Eq. (TTSl) we obtain that Ts s'{t,t') := Ts(t)Ts'{t) is always an 
eigenvalue of Ts^s'(i,iO := Ts_s/(^C, A^",A)';) with a product-singlet eigenvector |4'o) = IV'o)® IV'o)- One can further 
show that it is always a dominating and non-degenerate eigenvalue by demonstrating that — Ts{t)Ts' (t)l 

is a negative definite operator on Ha \ C|'I'o) (see Sects. A,B of [21| for details. We note though that we managed 
to rigorously prove negativity only for s < sq = 3/2 and further succeeded to confirm it with exact numerical 
computations up to much larger maximal auxiliary spin sq, while for any s it formally remains a conjecture). Denoting 
by Ts,s' (/r, A, /r', A') the continuation of the dominating eigenvalue in the proximity of the domain 'D~ x T>^, and using 
Hellmann-Feynman theorem to evaluate its first derivatives d^-Ts^s'in , , X^,) = ,X^)Tg,{fi^,,X^,), 


dxX 


/{fXf , Aj ,fi^,X^) = > Aj )9;^+r°,(/r/;, A/i), the HSK can be computed as 


Ks,s'{tH') = n[Ts{t)Ts>{t')] {ts,s'{h At > > A))) - t° (^i* At As'AtAt)) +0{e '>'”). 


(18) 


Remarkably, A term exactly cancels, while the finite- 
size corrections are exponentially small in the gap 7 = 
log |Ts,s'(t, t')/''"^! > 0 to subleading eigenvalue of r' of 
Ts_s'(t,tO- We shall later derive an explicit expression 
for HSK. 

What remains to be shown is that Xs{t) have well 
defined expansions in terms of local operators in TL 
n —> 00 . For any fc-local basis operator < 7 ^,. := 

ct/V 2 ^ • tr))'', ai^k ^ 0 , we write the component of 

O as [Ts(t)]~^dx+iaf:j,,Ts{nt)TsAt))- For treating 
n —>■ 00 asymptotics we substitute [L+°(t)/Ts(t)]”“^ = 
|' 0 o)(' 0 ol + 0{e~^^) and take into account the fact that 
the A—derivative should always hit the last, A:—th factor, 
producing SaLs = Sa^, otherwise the whole term would 
vanish due to the Eqs. (na. Thus we find a compact ma¬ 
trix product formula for the components (with the k = 1 
component vanishing) 

(a^„Xg(f)) = (V^„JX“=...X“-HV'a.)+0(e-^"), (19) 

where X“ := L+“(t)/rg(t), li/a) := V2is“JV'o)/T-s(A)- 
The HS norm of Xs{t) projected onto £ sites, in the limit 
n — £ ^ 00 , can be written analogously to Eq. (na 

e 

+ = ( 20 ) 

k—2 a 

(('>^ 0 1 At , At", , A+)] Vo) 

-(V'ol [AiArAr)YAo)Ao\ [L°(/r(^,A)^)]^|V'o)^- 

thus resulting in expression oc cf. Eq. (na, without any 
finite-size (^-dependent) corrections as |d>o) is an exact 
eigenstate. We have thus shown that the expansion 




k—2 a 


x—0 


( 21 ) 


is complete in the HS norm. Q.E.D. 


Eqs. (I20I21|) have two useful implications: (i) As 
the state |4'o) is a spin singlet (in 4—spin auxiliary 
space) the only relevant part of the SU{2) invari¬ 
ant TM Ts,s'At') = EaLj“(t) 0 L+“(f'), is the 
(2 J -|- l)-dimensional block, J = min{s, s'}, constitut¬ 
ing the spin singlet subspace of Ha, where it can be 
written explicitly as a tridiagonal matrix (see Sect A 
of 2l|). (ii) The HSK can be compactly written 
in terms of the resolvent of the TM, similarly as in 
[Hi, namely ^')] ^1^)) 

where Ts,s'{t,t') = Tg,g-(t,t')/[Tg(t)Tg/(t')] and |4') = 
Sae{x.y.z} IV'a) ® AA^ o.g. via solving a system of 2J 
linear equations 


Kg,g,(f,t') = n(T|$), (1-Tg,g.(i,0)l^) = 1^)- (22) 


By deriving the explicit form of matrix elements of 
Tg^s/(t,f') and solving Eq. (l22l) . we can encode the HSK 
explicitly in terms of a superposition of Cauch y-L orentz 
distributions (assuming s < s') [see Sect. A of |2l| ] 


Ks,s’{t,t') = 


Ks,s’ jt - t') 

Ts{t)Ts'{t') ’ 


(23) 




l{l + 2(s'-s))(2s + 1 - l)(2s' + 1 + 1) 


1=1 


(2s-bl)(2s'-f 1) 


" Cs' —s+Z ? 


where Cs{t) := 


s^ + A' 


Note that HSK is symmetric Ks,s'{t,t') = Ks',s{t',t) and 
strictly positive Ks,s'{t,t') > 0, Vs,s',t,f'. 

We would like to remind the reader that in TL n ^ 00 
the s = ^ family Xi (t) is equivalent to the family of local 
charges Qfc, as follows from Eqs. (17151) 


fk 

V(^) = E Qk+2 = d^Xi{t)\t=o- (24) 

fc=0 

Eq. (IT 9 )) thus generates also a handy explicit matrix prod¬ 
uct representations of the standard local conservation 
laws Qk or their densities qk- 
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Proof of linear independence.- Let us first show that 
Xi{t) are linearly independent from Xi(t), i.e., from Qk- 
We define an operator 


= / dt'Mt')xxt'), 


(25) 


where the function ft{t') is determined by minimizing the 
HS norm ||Xi(t)|||g, i.e. by the variation 

^^(Xi(t),Xi(t))=0, (26) 

resulting in the Fredholm equation of the first kind 

/ OO 

-OO 

Using the fact that the kernels (051) are related to Cauchy- 
Lorentz distributions Cs{t) up to trivial rescalings, we 
make an ansatz ft{t') = {Ti{t')/Ti{t))(p{t — t') which 

maps (EH) to a linear convolution equation |ci * v? = 
|c3, which, using the well-known convolution identity 

Cs * Cs' = TTCs+s', results in = ^ci, or 


Mt') = 


9^((3/2)2+t2)((i/2)2 + (i_i/)2)- 


(28) 


The conclusion of this analysis is that a family Xi{t) is 
(a) quasilocal (see Sect. D of 2l| for a numerical exam¬ 
ple) as its HSK, computed via Eqs. (I23I25I28L is exten¬ 
sive {Xl{t),Xl{t')) = ^,(t)",(tq (|ci(t-tO- ^C 2 (t-f')) 
and (b) is orthogonal to (and hence linearly indepen¬ 
dent from) all known local operators contained in the 
s = 1/2 family Xi/2{t), namely we have {Xi{t),Qk) = 
{Xi{t),Xi{t')) = 0, for all More generally, 

one can orthogonalize Xs(t) for higher s to all previous 
Xgr{t') for s' < s, by making an ansatz Xg{t) = Xg{t) — 

JZo with 

explicit expressions for bounded integrable functions 
fl^g-i(t'), fl g_i{t')- These families are HS orthogonal 

for different auxiliary spins, namely (Xs(t), (t')) = 0 

for s ^ s', while {Xg{t), Xg{t)) > 0, i.e. Xg{t) ^ 0, for 
all s,t (for details see Sect. C of [HI). This implies that 
Xg(f) are linearly independent from all previous Xgiif'), 
for s' < s, and in particular from Xi{t') or Qk- Q.E.D. 

Discussion.- We have proposed a direct extension of 
local conserved operators derived from the logarithm of 
the fundamental TM ia [Hj | to higher spin auxiliary 
spaces. We have proved that in such a case, the resulting 
operators are quasilocal. An interesting side-result of our 
statement is an asymptotic (thermodynamic), n —^ oo, 
inversion formula (^j T“^(^-|-it) ~ [Tg(t)]~^Tg{—^+it), 
valid for any s € which can be proven by im¬ 

plementing our matrix product formula m together 


with the gap statements (Sect. B of [2l|) to show that 
Tg{p.^)Tg{Xf) ~ rs(<)l. Our quasilocal operators Xg(t) 
© can thus be understood as logarithmic derivatives of 
Tg (Xf ). In TL n — > OO they become Hermitian opera¬ 
tors for any t G M.. For s = i, the Taylor expansion 
coefficients in t turn out to be local operators, while for 
s > ^, they remain non-local but quasilocal. One could 
thus equivalently work with a discrete series of quasilocal 
operators Qg^k +2 = {l/k\)d^Xg{t)\t=o, s € \Z+,k G Z+, 
rather than with a series of continuous families Xg (t). As 
a double index suggests, the number of relevant quasi¬ 
local charges in a large finite system may grow as 
rather than n as in the ultralocal case, although this ques¬ 
tion cannot be made precise with the results at hand. 


Our results promise a number of timely applications 
and generalizations. The new quasilocal families should 
be included in order to correctly describe fc —>■ 0,a; —>■ 0 
limit of dynamical structure factors and general Drude 
weights at finite temperatures [23 - E6| . or GGE in quan¬ 
tum quench protocols Q. For computing stationary ex¬ 
pectations of local observables after a quench from a non- 
thermal initial state, such as e.g. the Neel state |N), one 
can readily demonstrate extensivity (N|As(t)|N) oc n by 
extracting the leading eigenvalue of an associated transfer 
matrix, essentiall y p roceeding along the lines of calcula¬ 
tion done in Ref. jl^ for the fundamental (s = 1/2) TM. 
Appropriate g—deformations of the concepts developed 
in this Letter should provide additional quasilocal opera¬ 
tor families for the anisotropic Heisenberg model {XXZ 
chain). Extensions to SU{N) symmetric integrable spin 
chains seem straightforward, whereas a generalization to 
continuous quantum integrable systems and field theo¬ 
ries (such as Lieb-Liniger or sine-Gordon models) should 
be a challenge for the future. We close by stressing 
an important point of distinction with respect to spin- 
reversal symmetry brea king quasilocal conserved opera¬ 
tors in XXZ model 14l-ll7l|. Quasilocality, as abstractly 
formulated here, requires a finite-dimensional (but non¬ 
fundamental) representation of a quantum TM, and a 
factorizability condition for the leading eigenvalue of the 
associated auxiliary TM. This can happen, either for ir¬ 
reducible unitary representations of the symmetry group, 
but will result in operators which are always even under 
spin reversal, as is the case here; or due to the root-of- 
unity (commensurability) condition for the anisotropy, 
where highest-weight type non-unitary representations 
become reducible to finite dimensional ones, such as in 
the XXZ model. 


E. 1. thanks J. de Nardis and J.-S. Caux while T. P. 
thanks M. Mierzejewski and P. Prelovsek for valuable dis¬ 
cussions on related topics. M. M. and T. P. acknowledge 
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A: EXPLICIT COMPUTATION OF THE 
HILBERT-SCHMIDT KERNEL 

Let us define the following operator over the 4—spin 
auxiliary space 'Ha. = (2)fc=i 

■= - Ts{t)Ts'{t')l (SI) 

where has been defined in the main text. For 

s = s' and t = t', the operator Ts,s{t,t) is the TM for 
computation of the HS norm, so its maximal eigenvalue 
has to be positive as the HS norm is always positive. 
To complete the proof of quasilocality, it amounts to 
show that ¥s^s{t,t) is negative definite on 'Ha \ Cl^ho) 
where [I'o) = IV’o)® IV’o) is the leading eigenvector, which 
implies that [Ts(t)]^ is the maximal and non-degenerate 
eigenvalue with a finite gap to the sub-leading eigenvalue. 


This will be shown by explicitly constructing the operator 
(ISII) in a convenient basis (see section B for completing 
the proof of the gap), which in turn is needed for explic¬ 
itly computing the HSK Ks^s>{t,t') = Xs'{t')) via 

Eq. (12^ of the main text. 

Using the expansion 

a+b<2 

Fs..'(LO= E (S2) 

a.be{0,l,2} 

and suppressing dependence s, s' in ^ for compactness 
of notation, the non-vanishing matrix-components can be 
readily represented in terms of SU (2) invariant tensors 
IFa.h by employing a shorthand notation for elementary 
SU{2) symmetric operators over 'Ha- |*, jl := Sa^ • Sa^ 
and li,j, fc] := i(sa, x SaJ • Sa^, namely: 


Fo,o — — ^^(s + 1) — + 1) ~ s{s + l)s'(s' + 1) + |1,2]|3,4] — [1,3]|2,4] -|- |1,4]|2, 3] 

+ i ([1, 2, 31 - [1, 2,41 - [1, 3,41 + [2, 3,41) + i (|1,41 + [2,31 - [1, 21 - [1, 31 - [2,41 - [3,41), (S3) 


Fi,o = I ([1,31 + [2, 31 - [1,41 - |2,41) + i ([1, 3,41 + [2, 3,41), (S4) 

Fo,i = I ([2,31 + [2,41 - [1,31 - [1,41) + i ([1,2,31 + [1,2,41), (S5) 

Fi,i=-[1,31-I1,4]-|2,3]-[2,41, (S6) 

Fo,2 = -s(s + 1)-[1,21, (S7) 

F2,o = -s'(s' + 1)-I3,4]. (S8) 


Reduction to the singlet snbspace 


By virtue of SU (2) invariance of {t, t') the computation of HSK can be facilitated (assuming s < s' throughout 
this section without loss of generality) in the invariant 2s -I- 1 dimensional subspace formed by singlet eigenstates, 
'Ha U Vo = lsp{|j); j = 0,1,..., 2s}, which can be expanded explicitly in terms of computational basis |mim 2 TOlm 2 ) = 
|«T-i)ai|«^ 2 )a 2 l"i})a 3 |TO 2 )a 4 with help of Wigner 3j-symbols 


\j) = E E 

M=-j 


mi,7712 —— s 


/ S S 
\mi 7712 



^ 177717712777} m 2 ), 


(S9) 


where only the extremal singlet state factorizes |0) = 
I'l'o) = |'0o)aia2l'0o)a3a4- Let US denote the restriction to 
the singlet subspace as Fg°],(t,t') := Fs_g<(t,t')|vo. 

The proof is based on showing the following elemen¬ 
tary statements: (i) F^°], (t, t') is a quadratic form in the 
difference variable t := t — t , with identically vanishing 


linear terms F^^g = Fq°| = 0, i.e. 

fJ/ (t, t') ^ F®, (r) = D + F^^ (SIO) 
where D = F^^^ = F^^^ = 

(ii) is a real symmetric tridiagonal matrix in 

orthonormal singlet basis {|j)}. 

Both statements follow from demonstrating, by em¬ 
ploying certain elementary symmetry-based reductions. 
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that all matrices are expressible solely in terms of 
knowing only five types of matrix elements, namely 

01Il,3]|j), (Sll) 

(j|[l,3]|j + l), (S12) 

0-|[l,3][2,4l|j + l), (S13) 

(i|[l,2,3l|j + l), (j|[2,3,4l|j + l), (S14) 

for all \j) G Vo- The reductions of above expressions 
can be carried out after noticing simple transformation 
properties of |j) with respect to permutation operators 
Pi2 and P34 {Pij € End('Ha) swaps "Hai and "HaJ 

P 12 IJ) = P34|j) = (-i)^'+^'-'|j). (S15) 

We find straightforward implications 

0111,411/) = (j|l2,3]|/), (S16) 

0111,311/) = 0I[2,4]|/) (S17) 

0|[l,2,3]|/) =-0111,2,411/), (SIS) 

0l[l,3,4]|/) =-0112,3,411/). (S19) 

I 


for any //, where, in addition, a sign reversal under 
odd permutation of factors has been used for the triple¬ 
product terms. Notably, these identification are enough 
to establish vanishing of linear terms, = Fg°j = 0. 
Furthermore, with an aid of Casimir invariants 

[1,211/ = (-sO + 1) + 5 /j + 1)) I/, (S 20 ) 

[3,4l|j) = (-s'(s' + 1) + IjO + 1)) |J), (S21) 

we have Oll^a^li) = Oll^^^oli) = “l/j + !)■ With assis¬ 
tance of symbolic algebra in Mathematica, using explicit 
form of singlet eigenstates (IS^ for general s, s' and / we 
have been able to obtain 0|®'i°j|j) = j{j + 1), whereas 
vanishing of the upper-diagonal follows from Eqs. (IS16I) - 
(IS19I) . The T^-dependence (i) is at the end given by ele¬ 
ments (j|F[,°^|j) = 2(j||l,3l|j) = -i/j -k 1). 

Turning attention to the remaining (constant) term 
Fq^^q we first note that all triple-product terms vanish 
on Vo after resorting to explicit evaluation of matrix ele¬ 
ments from Eq. ()S14|) 


(j|[l,2,3l|j + l) 


-(j|[2,3,4]|i + l) 


U + 1)^ / (2s - j)(2P - /(2(s + 1) + j)(2(P + 1) + j) 

4 V (2j + l)(2j + 3) 


(S22) 


in conjunction with Eqs. (|S18IS19p . Notably, diagonal matrix elements of all triple-product terms vanish in accordance 
with their skew-symmetric nature. Furthermore, after repeating symmetry arguments based on Eq. (IS15P we arrive 
at the following explicit expressions: (a) for the diagonal matrix elements Uj = (j|Fg°/|/ 

= “i/i + 1) - s('S + l)s'(s' + 1) + + 1) - s(s + 1)) {^j{j + 1) - s'(s' + 1)) , (S23) 

and (b) for the first off-diagonal elements bj = (j|FQ°/|j -I- 1) = {j + l|FQ°g|j) 


/ = -(jII1,3]+2[1,3][2,4]|j + 1) = - 


j{j + l){j + 2) (2s - j)(2s' - j)(2(s + 1) + j)(2(s' + 1) + j) 


(2j + 1)(2 j + 3) 


(S24) 


We still owe the reader a brief remark in order to close (ii): by inspecting the structure of singlet eigenstates we 
can easily see that SU(2) invariant appearing in the expansion of Ts^s/(t,t') can raise/lower any magnetic quantum 
numbers at most by two, thus (t, t') is a banded matrix which cannot have non-vanishing matrix elements beyond 
the second upper/lower diagonals. At last, due to symmetry cancellations of the second upper diagonal between 
|1,3J|2,4] and |1,4]|2,3] projected onto singlets |/, we finally remain with a strictly tridiagonal form. Finally, 
F° is a symmetric matrix on Vq, i.e. (j|F° = (/|F° as consequence of mutual cancellation 

of skew-symmetric triple-product terms |i,/ fcj. 

It is useful to note that the matrix of the linear system given by Eq. (|2^ in fact coincides with F^°], (r), 

(1 -f,,,/t,t')) = (S25) 

Furthermore, note that the state |0) does not couple to the rest of 2s-dimensional singlet space Vq = lsp{|j),j = 
1,..., 2s}, since oq = &o = 0- Thus, by solving the 2s-dimensional tridiagonal linear system 

Fi°],(r)|S) = |1), (S26) 

we can readily obtain few explicit results for HSKs, 

with «:s,s/(t) =-^s(s -I- l)s'(s' -I- l)(l|s), (S27) 
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where the prefactor in front of Ks,s'(t) comes from expressing the state j'h) (see main text), expressed as |d') = 
(2/-\/3)-ys(s'+T)s^(sM-T)|l). More explicitly, and using a suitable gauge transformation (redefinition of bra-ket 
basis in Vq in order to remove the square-roots from off-diagonal matrix elements) and homogenizing the system 
(1^ . one can encode the HSK as 

Ks,s'(t) = -8s(s -I- l)s'(s' -I- 1) — , (S28) 

Xo 

where Xj satisfies a 3-point recurrence relation 


(j + l)(2s-b 1 -f j)(2s' -b l+j)xj-i + ji2s - j){2s' - j)xj+i + (2j -b l)iz-j{j + l))xj = 0, (S29) 


where z = 2 (t^ -b (s -b |)^ -b {s' + 5 )^), which can be solved with a direct backward iteration by choosing the initial 
conditions X 2 s-i-i = 0 , X 2 s = 1 - 

The solution, which is easily obtained explicitly for essentially arbitrary large s, can be neatly written in terms 
of a superposition of Cauchy distributions 


s{t) = 




(S30) 


and for few smallest auxiliary spins reads 


Ki.i 


2s(s-b 1) 
2s-bl 
8 20 

15 




-C2, 


K3 3 = —Cl -b 3c2 
2>2 16 


24 


84 


5 

= xC3 

3 = 


16^^’ 




63 

16^ 
144 


3c5 

2 

9 


12 56 28 20 

«;i.2 - yC2 -b —C3, - yC5 -b yC7, 


C 3 Q — ^ C 3 

2 .^ 52 

144 


■C3 


25 


■C4, 


21 


16 




21 24 

52 52 ’ 

28 24 36 20 

’5 155^52^52^32 


45 


80 

«^l,3 = ^C3 


36 


-C4, 


(S31) 


Moreover, a simple form for the superposition coefficients for all small s,s' lead us to conjecture that they can be 
written as low-order rational expressions. Indeed we found a remarkably simple closed form expression 


2s 


Ks,s’{t) = 


1=1 


l{l + 2{s'-s)){2s -b 1 - l){2s' + 1 + 1) 
(2s-bl)(2s'-bl) 


Cs'—s+liT), 


(S32) 


which reproduces the solution of Eqs. (IS28IS29I ) for any finite s, s', while we leave its rigorous derivation for the future. 

For a curiosity, we may write another closed form expression of HSKs for general s, s' and r, written in terms of a 
complex continuation of harmonic numbers known as the digamma function '4’{z)., 


m = iiogr(.) = E (i - Tzhi) - 

where 7 em is the Euler-Mascheroni constant, as 

,, s(s'(s' + l) + (s + l)2 + r2) ((,'-,)2 + + + 1 ) 2 +^ 2 ) ^ ^ 

2s+ 1 2(2s + l)(2s' + l) 


(S33) 


(S34) 


making use of a compact notation := ip{s -b s' -b 1 ± ir) — i/>(s' — s -b 1 ± ir). Remarkably, the second term gives 
rise to Cauchy distributions in (|S32I) via recurrence formula 


V'(z -b iV) - ip{z) = 


N gN, 


(S35) 


yielding 


^^s' + V's.s' = X! 


k^O 


1 


A: + (s' — s H- 1 + ir) A: + (s' — s + 1 — ir) 


— 2 Cs'('^) • 


k^l 


(S36) 
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B: FINITENESS OF THE GAP FOR THE AUXILIARY TRANSFER MATRIX 

Extensive ~ n scaling of a general HSK can be attributed to the finite spectral gap with respect to the leading- 
modulus eigenvalue of the entire Ha,- Thanks to Cauchy-Schwartz inequality 

Ks,s' (t, t') < (S37) 

it is sufficient to focus on s' = s and t' = t case only (pertaining to HS norm of As(t)), where we expand 

(S38) 

with matrix-valued coefficients reading 

A = -2s(s + 1) - ([1,21 + [1,31 + [1,41 + [2,31 + [2,31 + [3,4]) = -icW. (S39) 


1 = -i([l,4l - |2,31 - |1,2,3] - |1,2,4] - [1, 
C = -is(s + 1) - s2(s+ 1)2 + |1,21[3,4] - [1 
+i([l, 4] + [2, 31 - [1, 21 - [1, 31 - [2,4] - |3, 


These are just specialization of expressions given by 
Eas. lfSlVISS)) . The operator := (Sa^ -I- -|- Sag -I- 

denotes the four-fold s-spin Casimir invariant with 
eigenvalues s(s -I- 1). Note that the auxiliary operator 
denoted by C, such as in Eqs. (IS38IS4ip . should not be 
confused with a set of complex numbers. 

Denoting temporarily Ts^s{t,t) —>■ T{t) we note a re¬ 
markable commutativity property, 

[T(t),T(t')]=0, (S42) 

which is a direct consequence of Yang-Baxter equation. 
Specifically, considering a periodic chain of four spins s 
the auxiliary TM T{t) becomes the standard commuting 
quantum TM for the physical spin 1/2 now playing the 
role of auxiliary spin. This implies commutativity of all 
operator valued coefficients, 

[A, 1] = [A, C] = [B, C] = 0. (S43) 

In order to prove strict negativity of ¥s,s{t,t) on Ha \ 
C|4'o) it is enough to show that a quadratic (int) equa¬ 
tion ($|Fs^s(t,t)|$) = 0 does not have a solution, for 
any 1$) other than ITq). Due to (IS43I) this amounts to 
demonstrate that a matrix-valued discriminant 

A := l2 - 4AC, (S44) 

has only non-positive eigenvalues, while for any eigenvec¬ 
tor |$o) of A corresponding to zero eigenvalue, it must 
hold that (<I’o|Fs^s(t,t)|$o) < 0. 

Indeed, the entire singlet subspace Vo has the latter 
property, since we have A|vo = 0 due to A|vo = B|vo = 0. 
The negativity of Fs°i (t = 0) = Fg°Q on Vq follows from 
diagonal dominance of the tridiagonal matrix 

-aj>\bj\ + \bj_il j > 1, (S45) 


3,4]-[2,3,41), (S40) 

,3][2,41 + [1,4|[2,31 

+ i(|l, 2, 31 - |1, 2,41 - [1, 3,41 + [2, 3,4]). (S41) 


based on explicit form of matrix elements (IS23IS24I) . 

Clearly, for large enough t (Casimir) coefficient A 
starts to dominate and therefore (non-singlet) eigenstates 
belonging to any higher spin multiplet necessarily be¬ 
come sub-leading and the spectral gap 7 > 0 is always 
due to the largest (smallest in modulus) (singlet) eigen¬ 
value of Fg°g. For a generic t € R on the other hand it 
might happen that the gap 7 is determined by eigenvec¬ 
tors outside of Vq. At the moment we have only been 
able to rigorously confirm our statement for s € {i, 1, |} 
by analytically diagonalizing the operator A projected 
onto highest-weight total spin S' > 0 subspaces of Ha 
{SU{2) descendants only contribute to degeneracies), or 
some larger s, by extensive numerical checks. 


C: FREDHOLM-GRAM-SCHMIDT 
ORTHOGONALIZATION FOR HIGHER 
AUXILIARY SPINS 

Using the appealing explicit form of HSK (1^ . derived 
in Sect. A, we here outline a general scheme of orthogo- 
nalization of Xs{t) to Xs'{t') for all s' < s, t' € R. We 
denote such orthogonalized quasilocal conserved opera¬ 
tors as Xs{t). Picking a set of suitable functions // j,/(t'), 
for s' G iZ+ < s, we seek for an operator 

^ poo 

Xa{t) = Xs{t) - ^ / dt'flAt')XAt'), (S46) 

s' 

which minimizes the HS norm ||Xs(t)||yg, i.e., 

= 0, s' < s. (S47) 
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This yields a coupled linear system of (2s — 1) x (2s — 1) 
Fredholm equations of the first kind 


to (2s — 1) X (2s — 1) system of linear equations with 
coefficients that are polynomials in variable z 


^ poo 

{f, {t") = 

J-oo 

Vt', s' < s. (S48) 

If the unknown functions /* gi(t') are sought for in terms 
of the following difference ansatz 

(S49) 

then, noting that the HSK (23 18271) also obeys a scaled 
difference form 

Ks',s"{t' - t") = n~^Ts’{t')Ts”it”)Ks’,s"{t',t''), (S50) 

the Fredholm system (I848|) becomes equivalent to a linear 
convolution system 

s” <S 

Ks',s" * V’s^.s = Ks',s, s' < S, (S51) 

s" 

where {ip * ip'){t) := dt'ip{t')if'{t - t'). 

For s = 1, this yields a single condition 


s" <S 

Crs',e"{z)gs',s{z) = Gs',e{z), 
s" 


(856) 


where (now assuming s < s' without loss of generality) 
Gs,s’{z) ■■= A(ks^s/)(z) = Gs’,s{z) = (857) 

((( + 2(s'-s))(2s + 1 - l){2s' + 1 + () 2 (s'-s +0 
^ (2s + 1)(2s' + 1) 


Elementary algebra yields a solution which is nonvan¬ 
ishing only for the last two components s' = s — 1 and 
s' = s — i (note that here s > 1): 


gs',s{z) = 0, for s' < s - 1, 

1 r / -1 

gs-iA^) = 




(858) 

s(25—1) 2 

s{2s-l)-l^ ’ 


9s-ls{z) 

where 


, 1 _ 

(2s-l)(2s-fl) s(2s-fl) V 

s(2 s-bl) 


(5_1)(2s-1)- 



(859) 


-Cl * =-C3, (852) 

with a unique solution, equivalent to (1281) . 

= (853) 

which follows directly from an elementary addition iden¬ 
tity for the Cauchy distributions 

Cs * Cs'= Trcs+s>. (854) 

For s > 1 the system (18481) becomes nontrivial. Then, 
it turns advantageous to construct a linear isomorphism 
A : ip ^ g between the convolution ring of integrable 
functions (or distributions) ip{t) with operations (+, *) 
spanned by {cs, s G ^2)+}, and the ring of functions g{z) 
of a formal variable z with operations (+,•)) where • is 
the usual pointwise multiplicaiton, analytic on the unit 
disc around the origin z = 0. The map A and its inverse 
A~^ are defined uniquely by: 

A{cs) = ttz^% A{ip*ip') = A{ip)A{ip'), (855) 

A-yz'') = -Ck, A-'^{gg') = A’^p) * A'^^'), 

TT 2 

and the linearity. Note that a constant function in the 
image of A corresponds to a Dirac distribution co(t) = 
d(5(t), which, however, never appears in our calculation. 

Using the following notation for the unknown functions 
gs',s '■= A{ipsps), the Fredholm system (IS48I) is A-mapped 


Note that the convergence radius is always larger 
than 1, guaranteeing analyticity inside the unit disc. Ex¬ 
panding the geometric series and transforming back with 
A“^ (18551) . we obtain explicit results for the two nonva¬ 
nishing functions 


/L-i(i') = 


( s{2s-l) 

oo \ 


rs-i(t') ( s(2 s-l) 

7rrs(t) 


(s-1) 


1=1 


2{2sAs_i{t') ( 


1 


s(2 s-|-l) 


(S60) 


1=0 


which complete the explicit construction of Xs{t) (18461) . 
We note that the exponentially convergent sums above, 
Eqs. (18601) . allow closed form expressions in terms of the 
Hypergeometric function 2 A 1 , or the incomplete Beta 
function, of argument 1/^s and with complex parame¬ 
ters. 

It may be of interest also to consider HS-norms and H8- 
kernels defined with respect to orthogonalized quasilocal 
operators 


Ks{t,t') 


{Xs{t),Xs{t')) 


Ks{t - t') 

71 - 

Ts{t)Ts'{t')' 


(861) 
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For example, showing that Ks{t,t) = ||^s(i)|lHS > 0 is 
a final step of the proof that the Xs{t) are linearly inde¬ 
pendent for different s. In the opposite case, specifically 
if, for some s, Xe{t) would be expressible as a linear 
combination of Xs'(t'), for s' < s, then one would have 
Xs{t) = 0, and hence Ks{t,t) = 0. 

Clearly, designating Gg = A(ks), we find 


s' ,s" <s 

Gg{z) = Gg^s{z) - ^ 9s',s{z)Gg>^g»{z)gs",s{z) 

s' ,s" 
s' <s 

= G,,,(z) - ^ gs'A^)Gs'Az) (S62) 


(2s)^7r 


(s-l)2(2s-l)(2s-bl) 


s(s—l)(2s —1) 

27+1 


2" + 1 - 1 - — 


and transforming back 

Ks{t) = 


{2sr 


(s-1 )2(2 s-1 )(2 s + 1 ) 

s(s— l)(2s— 1) 

-ciK'i )- 

1=1 


2s + l 


A) 


(S63) 


Specifically, noting that Cs'(O) = I/s': 

(2sA rs(s—l)(2s —1) 


c(0) = 


7r(s-l)2(2s-l)(2s + l) L 2s + l 
4s3 - 2s + 1 


log 


s(2s - l)(2s+ 1)J’ 


(S64) 


which satisfies k(0 ) > 0 for any s > 1, and hence 
||Xs(t)||Hs = Ks{t,t) = nK(0)/[Ts(t)]^ > 0. Note that 
the case s = 1 has been treated separately before. 


D: NUMERICAL EXAMPLE 


thogonal to all the local ones, namely Xs=i{t = 0). Ma¬ 
trix product formula (IT^ inserted to Eq. (1251) with (l28l) 
yields all the local terms in the infinite size limit n —>• 00 , 
say up to support size £ < 4 




n—1 


x—0 


155^ ^ 

* f^x-t-2 4” 252 


16 _ 53 _ ^ 

+ -^[Crx-<^x+l){0'x+2-0'x+3) — -^Ax-0'x+2)Ax+1-<^x+3) 


^^Ax'^x+3')Ax+1'^x+2')^ + h.O.t. 


(S65) 


We note that this qualitatively agrees with the optimal 
quasilocal conserved operator Q' which has been con¬ 
structed approximately by a systematic numerical pro¬ 
cedure in Ref. [1]. Small quantitative deviations in the 
coefficients (note that [ 1 ] used spin- 1/2 operators rather 
than Pauli matrices which attributes a relative factor 
of 4 in quartic/quadratic terms) can be explained by 
the fact that the operator (IS65I) . being just one mem¬ 
ber of the s = 1 family Xi{t), is not optimized with 
respect to a relative weight within a finite support i, 
Xe{X) = lim„^ooEfc=iEal(f^fe:^)lV(2^,2r). On the 
other hand, the operator Q' of [1] is determined pre¬ 
cisely by maximizing Xe{X) within a given set of con¬ 
served X. For Xi(0) we obtain = 0.508,0.682,0.797, 
for i = 4,5, 6 , respectively, while for the optimal nu¬ 
merical Q' one has [1] Xg = 0.605,0.759,0.840. More¬ 
over, numerical inspection of relative weights gg{X) = 
Xg{X) — Xe-i{X) of a sequence of higher quasilocal op¬ 
erators, A's(O), s = 1 , |, 2 , indicates that for larger s the 
relative weights fj,i{Xs{0)) have clear maxima at larger 
supports £ ^ £*{s), while after that they decay exponen¬ 
tially fig ~ e~'^',£ > £*(s). 


Here we write out, explicitly, the leading terms of the 
simplest new quasilocal conserved operator that is or- 


[1] M. Mierzejewski, P. Prelovsek, and T. Prosen, Phys. 
Rev. Lett. 114, 140601 (2015). 


















